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Abstract. We generalize the method introduced in J. Phys. A: Math. Gen. 35, 
7255 (2002) based on the concept of thermodynamic equivalence and we transform a Fermi 
system of general density of states into a thermodynamically equivalent Bose system. This 
consists of mapping configurations of fermions from the original system onto configurations 
of bosons, the initial and final configurations having the same energy above the many-body 
ground state energy. In this way we obtain two systems of particles of different exclusion 
statistics, but which have the same entropies-and therefore identical canonical thermody- 
namic properties. This method enables one in general to calculate the system properties in 
either of the bosonic and fermionic representations. We check the method here in micro- 
scopic detail by calculating the equilibrium particle distributions in the two representations 
using the entropy maximization at fixed particle number and fixed "fermionic" and "bosonic" 
energies, respectively. Analytical calculations seem difficult to do, but we check the results 
numerically and we find them equal within the numerical accuracy. 
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1 INTRODUCTION 

The heat capacity of ideal systems of constant density of single particle states 
(DOS) and fixed number of particles is independent of the exclusion statistics 
of the constituent particles [HEIEIIIIEIEIEIE]- Moreover, integrating the 
ratio between the heat capacity, C, and the temperature, T, with respect to 
the temperature, we obtain the entropy, S = Jq {C/T')dT' . So S, expressed 
as a function of T, particle number N, and all the other extensive parameters 
of the system (let us denote them by X), is also independent of the exclusion 
statistics. Based on this property, systems of equal, constant, DOS have been 
called thermodynamically equivalent [6]. 

The thermodynamic equivalence relation can be extended from the set 
of systems of constant DOS to the set of all physical systems [8]. To define the 
thermodynamic equivalence from a more general perspective, let me split the 
internal enegy of the system as U = Ugs + Ub, where Ugs is the energy of the 
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system at zero temperature and Ub is the excitation or Bose energy. Then two 
systems of entropies Si{Ui, N,X.) and S2{U2, N,X.) are thermodynamically 
equivalent if Si{Ub + Ug.s.,i{N, X), N, X) = S2{Ub + U^.s.M^, X), A^, X). The 
thermodynamic equivalence splits the set of all physical system into equiv- 
alence classes-like the classes which include the ideal systems of the same, 
constant DOS [HE]. 

The equality of the entropies of two systems with the same Ub and {N, X) 
implies that the number of microstates with the same excitation energy is 
the same in the two systems for any values of Ub and extensive parameters 
{N, X). This was shown explicitly to be true for systems of the same, constant 
DOS, by mapping configurations of particles of arbitrary exclusion statistics 
onto configurations of bosons of the same Ub [S] . Vice- versa, the mapping of 
configurations of ideal bosons into configurations of ideal fermions puts into 
correspondence the bosons on the ground state (the Bose-Einstein condensate) 
with a degenerate subsystem of fermions at the bottom of the single particle 
spectrum, which was called the Fermi condensate [HI [9]. The configurations 
of fermions which contain the Fermi condensate have higher probability than 
the "standard" Fermi distribution, as it was shown in Ref. [S], contrary to 
what was generally belived. In interacting systems, the existence of the Fermi 
condensate may modify so dramatically the weights of different configurations 
of particles that it may lead to first order phase transitions [U [TO] . 

In this paper we generalize the method of Ref. [8] and we show how 
to map a Fermi system of general DOS to a thermodynamically equivalent 
Bose system. The technique may be extended (although not in a trivial way) 
to any transformation between Bose, Haldane, and Fermi systems and, to 
distinguish it from other transformations, it will be called exclusion statistics 
transformation (EST). 

In Section 12.11 we calculate the quasiparticle energies and the density of 
states of the Bose gas. If the Fermi DOS is not constant, the Bose quasiparticle 
energies depend on the population of the other quasiparticle levels; in other 
words, the bosons are interacting. This effect is discussed in Section [3 

By EST, the particle population of the energy levels in the Bose system 
is determined by the particle population in the Fermi system and vice-versa. 
On the other hand, once we set the rules for calculating the DOS and the 
interaction between the particles in the Bose gas, the particle populations in 
the two systems (Bose and Fermi) can be calculated independently, by entropy 
maximization. In this way it is shown that one can make calculations in either 
of the systems and transfer the results to the other one. 

Last section is reserved for the discussion of the results and conclusions. 
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Figure 1: The illustration of the basic idea of exclusion statistics transforma- 
tion (EST). 



2 THE TRANSFORMATION FROM FERMI TO 
BOSE SYSTEMS 

Concretely, at microscopic scale the EST from a Fermi to a Bose system may 
be imagined as follows: we picture the fermions as forming groups of "close 
packed" particles that occupy completely certain intervals along the single 
particle energy axis (no empty single particle states in the intervals) and we 
associate to each of these groups an energy level in the equivalent Bose system 
(see figure [T]). The groups are labeled with i = 0,1,..., starting from the 
bottom of the energy axis and going up. So, the group i of Ni fermions 
is associated to the level eB,j of Ni bosons, in such a way that eB,o = 0, 
eB,i ^ cbj if ^ < j; and the total energy, E-q^ of the new system is the same 
as the excitation energy of the original Fermi system-the excitation energy of 
the Fermi system being the difference between the energy of the system, E-p, 
and the lowest energy attainable at the given particle number, Ugs (in other 
words, C/gs is the energy of the system at zero K). 

To distinguish between the energies of individual configurations and the 
(average) internal energy of the system in given thermodynamic conditions, 
we denote by E-p the energy of a configuration of fermions (the Fermi energy 
will be denoted by ep) and by E-q{= Ep — Ugs) the energy of a configuration 
of bosons. Since for a given excitation energy, £'b, to each configuration of 
fermions it corresponds a configuration of bosons and vice-versa, then the en- 
tropy of the Fermi system, Sp, of internal energy U = Ugs + EB, is equal to the 
entropy of the Bose system of internal energy Ub = Eb [Sp{Ugs + Ub, N, . . .) = 
Sb(Ub, N, . . .)], i.e. the two systems are constructed to be thermodynamically 
equivalent. With these definitions we have = SSb/OEb = dSp/dEp 
and fiB/T = -dSB/dN = (/i - ep)/T, where /i/T = -dSp/dN and ep = 
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dUg,s./dN [8\. In figured] we see also how the Fermi condensate ^ becomes 
Bose-Einstein condensate in the Bose system [8l [11] . 

2.1 BOSONIC SINGLE PARTICLE ENERGIES 

Let's now calculate the bosonic energy levels. We shall use e and eb to denote 
the single particle energies in the Fermi and Bose systems, respectively. If m 
is the number of particles in the system, the particle number dependent Fermi 
energy is denoted by ep{m). Moreover, in both (Bose and Fermi) systems we 
number the single particle states form zero to infinity, starting from the lowest 
energy level and going up- wards; the order in which we number the degenerate 
levels is not important. At zero temperature, the fermions occupy all the 
states, from zero to the Fermi energy, €f{N). This is the ground state (g.s.) 
of the Fermi system and is put into correspondence with the ground state of 
the Bose system, ee.O) like in figure [2] (a). 

The correspondence between the Fermi and the Bose energy levels is 
done recursively, starting from the lowest energies, as it is depicted in Figs. [2] 
(b), (c), and (d). If the lowest hole in the Fermi system is created by lifting 
a particle from the level i, of energy €i{< ep{N)) to the level + 1, then in 
the Bose system i — 1 particles are left on the g.s. and N — i + 1 particles are 
lifted on the first excited level, eB,i [figure [2] (b)]. The excitation energy in 
both systems is the same, so we define 

_ €FiN + l)-€i 
£3,1 = 



N-i + l 



The recursion procedure continues upwards; let's say that we arrived at the 
level €j. We assume that we lifted k particles from the energy levels below €j, 
so we created already k quasiparticle energy levels in the Bose system. At least 
from ej, upwards, all the Fermi levels are occupied up to eF(A^ + A;) (i.e. the 
total number of particles, A^, plus the k particles that have been lifted above 
the Fermi energy) and all these particles stay on the highest occupied energy 
level in the Bose system, eB,fc- Let us now lift the n particles from the energy 
levels j + 1, . . . , j + n, onto A^ + A; + l,...,A^ + A: + n. This corresponds in the 
Bose system to n — 1 free states, say eB,k+i, ■ ■ ■ , ^B,k+n-i, and an n^^ state, 
eB,fc+n) with a population of A^ — {j — k) particles. The difference between the 
Bose levels which are successively occupied, eB,fc and eB,k+n, is 

6B,.+n-6B,. = N + k-j N + k-j 

(see figured]). 

The procedure continues upwards until we reach the highest energy par- 
ticles in the Fermi system. 
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Figure 2: The correspondence between the Fermi and Bose microscopic dis- 
tributions: (a) both systems are on the ground state; (b) four particles on 
the ground state, £3,01 and six particles on the first excited state, eB,i; (c) 
four particles on eB,0) three on eB,i, and three on ee.s (eB,2 is free); (d) final 
distribution, with four particles on eB,o, three on eB,i, zero on 63,2, two on 
eB,3, zero on 63,4, and one on £3,5- 



2.2 CONTINUOUS LIMIT 

To discuss the thermodynamic limit, we need to go from the discrete case to 
the (quasi) continuous limit. For this, we denote the fermionic and the bosonic 
DOSs by cr(e) and crB(eB), respectively. We assume that in the Fermi system 
A^o particles are condensed [9], so there are no holes in the energy interval 
from to eF(A''o). Above eF(-/Vo) we divide the single particle energy axis 
into microscopic intervals, numbered i = l,...,oo. Then the interval i of 
length 6ei, contains di = a{ei)Sei single particle states. The population of the 
interval 6ei is denoted as 6ni = f{ei)di and we assume that /(ej) < 1 for any 
i > 1. After EST, A'^o represents the number of particles on the Bose g.s., 
eB,0) and the intervals dti transform into the intervals de-Q.i- The number of 
bosonic states in each of these intervals is d-Q,i = di — 5ni + 1 ~ — 6ni (for 
large enough d-B,i)- Although it has no relevance in the continuous limit, note 
that because of an overlap between successive Bose intervals, the expression 
dB,i = di — 6ni, leads to a better counting of total number of bosonic states 
(the same argument is in [E]). Applying equation ([1]), 6eB,i is calculated by 
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exciting [1 — f{ei)]di fermions from the interval 5ei, which leads to 

From equation ([2]) we obtain the Bose DOS, 

cb Ub i = T — = ■■ U>) 

' S^B,^ eF[iV + 1 + E)=i{dj - 6nj)] - ^ ^ 

and the Bose population (the number of particles divided by the number of 
states) , 

o(eB,i) = -1 7— = -jzm 7 • ^) 

di - drii f ^[ei{eB,i)\ - 1 

Changing the summations into integrals, equation ^ becomes 

where e(eB) is the single particle energy in the original Fermi systems that 
corresponds to the Bose energy ee- 

To calculate e(eB)) note that the number of Fermi states in the small 
energy interval 5e which corresponds to the Bose energy interval Jee is equal 
to the number of particles in the interval, plus the number of bosonic energy 
levels: d{6e) = [b{eB) + l]crB(eB)'^eB- Integrating over all the Bose energy 
levels, we get the expression 

e(eB) = ep (a^o + H M^M^) + l]de| ■ (6) 



The total Bose energy, 

Eb= deBcrB(eB)eB&(eB) , (7) 
Jo 

is nothing but the excitation energy of the system, in agreement also with the 
definitions given in \E\. 



3 INTERACTION IN THE BOSE SYSTEM 

In the general case, the Bose system obtained after the EST is not an ideal 
system. A change in the population of Bose "single particle" states, say Sb{eB), 
changes all the single particle energies and, further, (Tb- To see this, let us 
insert ^ particles in the interval i of the Fermi system. By doing this, we also 
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Figure 3: Example of insertion of two particles (filled circles) on the fourth 
Bose energy level. 



expand the interval by ^ states, so that the number of Bose states in any of the 
intervals would not change (see figure [3]) . All the quantities that are affected 
by this insertion of quasiparticles will bear a superscript which refers to the 
energy interval (i), or directly to the energy levels [ej (Fermi) or eB,i (Bose)] 
where the insertion took place. Applying equation we calculate the new 
energy interval Se^^ y = e]^ . - ej^ 



5e 



B,fc 



eF[N + l + ^ + J2{dj - 6nj)] - ep 



or 



dk - 5nk 



No + J2dj 



iV + ^-A^o-ELoH 



if i > A; . 



If ^ <C 5ni, then 



d — 5nk 



de-p 



(8) 



(9) 



[N + l + Y^{dj-5nj)] 
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-^i^o+j:d,]\ , ifi<k (10) 

or 



x^[N + l + J2{d-6nk)], ifi>k. (11) 

If in the Bose system we add up the energies of all particles, we obtain 

= Eb + S,eB,i + Ej 

where Eb is given by equation ([7]) and E'l oc ^ is the interaction term. We 
may interpret ee as the quasiparticle energy, and write it as 

eB = CBf + X! ^(^Bf , e'Bf)dn{eB{) , (12) 

where eef is the energy of the "free particles", while v is the interaction po- 
tential. Note that v{eBi,e'B{) depends on n(eg£), for any value of eg£. 

3.1 CALCULATION OF THE INTERACTION POTENTIAL 

The Bose quasiparticle energies and the interaction potential v from equation 
p2p may be calculated by adding the microscopic energy intervals given by 
the equations ()10p and (jlip . Transforming the summations into integrals and 
using the equations ^ and ([6]) we obtain 



Jo 

JO 



N - No - S^"" h{e)aB{e)de 

[^ + l + /o"^B(e)d6]d6B 



eF[A^ + 1 + /o^« ^B(e)de] - ep {A^o + /o^ <yB{e)[b{e) + l]de} 



/o iV-iVo-/o^«6(e)cTB(e)de 
for eB,j > ee 



and 



,^B..,^ ^ ^ p. l-CTB(6B)i^[A^ + l + r aB(e)d6] 

= " ^ iV-iVo-/o-6(e)aB(e)de 
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ree dee {^[iV + 1 + /o^" Me)de] - ^[Np + /J^ a^ie)[b{e) + l]de]} 



eF[N + 1 + J^^ aB{e)de] - ep {iVo + /o^ <^B(e)[6(e) + l]de} 
^B. 1 - aB(eB)^[iV + 1 + J^^ aB(6)de] /-^b 



AT AT r^BU \ f \A deB + ^/ crB(eB) 



{^[iV + 1 + f^- <TB(e)d6] - ^[iVo + F ^b(6)[6(6) + l]d< 



^-^^o-/o^ 6(e)aB(e)de 



deB 



for eBA < CB- 



(14) 



In terms of e, and using the shorthand notation eo = ep^Np), the equations 
(fT3ll and ([HI) become 



e,.,^ , r(^«) a(6l)[l-/(6l)] 

eB(?) = eB-C/ dei ^ — --- — 

o-B(eB(ei)) 



N -Ni 







de2cr(e2)/(e2) 







(iV + /;M62a(e2)(l-/(e2))) 






iV + /;M62CT(e2)(l-/(e2)) 


- ei 



, for a > e(eB) 



(15) 



and 



^B (?) - CB - ? / dei — -— - — 

Jea o"B(eB(ei)) 



-1 





ep 


(iV + /;Me2a(62)(l-/(62))) 




ep 


iV + /;Me2a(e2)(l-/(e2)) 


- ei 



N-No- I de2a(e2)/(e2) 

£0 



f<:(<:B) 

+ e / deicT(6i)[l-/(ei)] 



ep (iV + d62 a(e2)(l - /(62)))] - a-\e,) 



iV-iVo-/;:de2a(e2)/(e2) 



, for ei < e(eB), (16) 



respectively. The terms proportional to ^ in equations (I13p . ()14p . (llSp . and 
represent the interaction term u(eB,eB,i) of equation (fT^ . 

3.2 BOSE DENSITY OF STATES 

The changes of the single particle energy levels produce changes of the Bose 
DOS, <TB(eB)- The new DOS is obtained by using eqs. ([TU]) and PT]) : 



<7B(eB,fc) 

dep 



d — 5nk 



o-B(eB,- 



B,fc 



1-C 
dep 



0"B(eB,A 



N-Np- E •=! <5n, 



(17) 



X < ^[N + l + Y^{d,-5n,)] - ^[Np + Y^d,] ^ I , if ^ < fe, 
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and 



f^B(eB,fc) = —r-i ~ o"B(eB,fc) • U + 



o"B(eB,fc) 



(18) 



dep 



[iV + 1 + ^(d-5nj)] , if i > fc. 



Writing equations (|T3|) . (jMj) . p5|) and p6j) in the format 



(19) 



and using equations (fTTI) and ([18]) we obtain 
See,/,. 



o-B(eB,fc) • 1 



r''''o-B(eB„ 



j dep dep 

jdiV [jV+l + /;B,fe.3(e^)d.^] div 

if eB,i < ^Bh, 



Ar_iVo-/o^«n(e'B)aB(ei3)dei3 

{iVo+/;BaB(e^)[fe{^^) + l]d.^}}) ' 

(20) 



and 



cTr(eB,;c+r«- 

X < 1 - crB(eB,A;) 



dep 
diV 



o-B(eB,fc) • 1 + 



Writing the integrals in terms of e, Eqs. (|20p and (|2ip become 



iV - iVo - /o"''= b{e'^)aB{e'^)de'^ 
] , if eB,i > eB,k . (21) 



0"b (eB,A:+? 



B,fc 



o-B(eB,fc) • 1 



C''''"o-B(eB,fc) 



Q^e^,' j^<|.deia(ei)/(ei) 
er (iV + dei a(ei)(l - /(ei))) - ^7-1(6^)^ 



= fTB(eB,fc)- |l-f«'^ 
if eB,i < eB,fc 





ep (iV + /;'=deia(6i)(l- 


-/(^i))) 






^iV + /,;'=deia(6i)(l 







(22) 



and 



■) = crB(eB,fe) • 1 + 



EB.i 



X O 



0"B(eB,fc) 



ep(iV + /;''deia(6i)(l-/(ei)) 



iV-iVo-/;'^d6ia(6i)/(ei) 
, if eB,i > EB.fe • (23) 
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4 FINDING THE EQUILIBRIUM DISTRIBUTION 



The Bose and the Fermi systems are thermodynamicahy equivalent by con- 
struction. Therefore, the equihbrium particle distribution in the Fermi gas, if 
transformed by ([!]), should give us the equilibrium particle distribution in the 
Bose gas and knowing this, we can calculate all the microscopic parameters 
of the Bose system and the quasiparticle-quasiparticle interactions like in sec- 
tions [2?2] and [STTl determining in this way all the properties of the gas. The 
Fermi equilibrium (i.e. most probable) particle distribution is the well known 
/(e) = [exp(/?e - /3/i) + 1]^^ (where /3 = (kBT)-^), and this gives us by dU, 
the bosonic population 

b{eB) = exp[/3;u - /?e(eB)] . (24) 

which can be readily used in the calculations. 

Nevertheless, we have still another method to calculate the equilibrium 
particle distribution in the Bose gas. In sections [2] and [3] we defined self- 
consistently the quasiparticle energy levels in the Bose gas and the quasiparticle- 
quasiparticle interaction, so we can use this as a starting point to build the 
Bose entropy, Sb, and canonical partition function, Zb- Then the equilibrium 
particle distribution in the Bose gas should be the one that maximizes Zb at 
constant and Ub- The Bose and the Fermi systems are thermodynamicahy 
equivalent, so the canonical partition function of the Fermi system, Zp, is 
identical to Zb by construction, and indeed, the distribution that maximizes 
one should correspond (by an EST transformation S]) to the distribution that 
maximized the other, unless the canonical and grandcanonical ensembles in at 
least one of the systems are not equivalent. 

To check this equivalence, let me calculate and maximize it at con- 
stant Ub and A'^. (Note that constant Ub and is equivalent to constant U 
and N.) So let again 5eB,k be the microscopic energy intervals along the ee 
axis, each containing dB,k states. The total number of configurations available 
in one interval is 

iSnk + dB,ky. , . 

WB,k = — ? — r; — r- ■ (25) 
onkldB,k'- 

Summing up the contributions of all the intervals and applying the Stirling 
approximation, logm! ~ mlog(m/e), we obtain the Bose entropy: 

5bAb=/ deBM^B){[l + b{€B)]\og[l + b{eB)]-b{eB)\og[b{eB)]}, (26) 
Jo 

To find the equilibrium configuration, we have to maximize 5b as a functional 
of b, at fixed Eb and N. For this we introduce the Lagrange multipliers (3b 
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and Pb/J-b, corresponding to Eb and N, respectively, and we solve 
6iS/kB - I3bEb + PbI^bN) 



6b 



0, (27) 



where, as shown in section [21 from thermodynamic equivalence arguments we 
have /3b = /3 and ^ub = /U — ep (see also Ref. [8]). The functional derivative 
6EB/6b{eB) is not simply eB- To calculate it we rewrite Eq. (fT2|) as 

ei' = eB+ de'^vieB,e'^)aB{e'B)5b{e'j,), (28) 
Jo 

where eB is the quasiparticle energy at equilibrium configuration and Sb(eB) 
is the deviation of the population from the equilibrium configuration, 6(eB)- 
Note that both, 6(eB) and 6b{eB) are given as functions of the energy before 
the insertion of particles, eB- The interaction term is calculated in Section [3. 11 
Using (p8]) we calculate the energy 

E^ = / de§ei'[aB{ei') + 6aB{ei')][b{eB) + db{eB)] (29) 
Jo 



roo r roo 

/ deB eB + / deB-u(eB,eB)crB(eB)56(eB) 
Jo L Jo 



M^bMcb) + Sb{eB)] , 



where we used the equation de^/deB = CB(eB)/[cB(eB ) + '^'^B(eB )] 
obvious notation 6(eB) = &[eB(eB )] ™ ^^'^^ line. From ([29]) . and keeping 
only the first order in 6b, we obtain 

POO 

E'^=Eb+ / deB eB^B(eB)56(eB) (30) 
Jo 

+ / deB/ deB6(eB)cJB(eB)w(eB,eB)c^B(eB)^^(eB) • 
Jo Jo 

The first integral in equation (|3U|). let's call it Eqp, is the typical contribution 
of the added quasiparticle energies to the total energy of the system, whereas 
the double integral (say E"^) is due to the collective change of quasiparticle 
energies and can not be eliminated by a redefinition of the ground state energy, 
as e.g. in the Fermi liquid theory. So E^ will produce an effective quasiparticle 
energy 

eB(eB) = eB+/ de'^b{e'Q)aB{e'Q)v{e'^, cb) = eB + eBieB) ■ (31) 
Jo 

Since 6S/6b\ej^ = A;BcrB(eB) log[l + 6"^(eB)] and 6N/6b\eg = o"B(eB), plugging 
(j3U|) into equation (|27p . one obtains the following equation for the equilibrium 
distribution: 

log[l + b-\eB)] - peBieB) + PfiB = 0. (32) 
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Equation ([32]) leads to the typical Bose quasiparticle level population, 

6(eB) = [exp(/3eB(eB) - /?/xb) - (33) 

If neither the Fermi, not the resultant Bose system is pathologic, then the 
expressions (|24p and (j33p should give the same result, eventually with some 
corrections which vanish in the thermodynamic limit. 



5 COMPARISON OF THE BOSE AND FERMI 
DESCRIPTIONS 

To make concrete calculations and check if Eqs. (j24p and (j33p give the same 
results, we assume that the Fermi DOS has the general form cj(e) = Ce^, where 
C and s are constants and we calculate the Bose energies, DOS, and chemical 
potential by using the Fermi distribution, /(e) = [exp(/3e — /?/x) + in the 
formulas of Section [2^21 Using these results and the expressions for v{e'-Q,e'B) 
calculated in Section [3l we can compare Eqs. (pll) and (p3]) . 



5.1 CONSTANT DENSITY OF STATES 



In Ref. [8] it is proven that the Bose and Fermi descriptions (obtained by 
EST) of systems with constant DOS are equivalent under canonical conditions. 
Therefore we start by using such systems to test the formalism introduced 
above. If the density of states is constant [a = C) then eF(ra) = n/C and by 
using ([6]) to rewrite of equation ([3]) as 



c- 



iVn 



1 + 

i=i 



the Bose DOS becomes [S] 

o'B(eB,j) = C = a . 
Moreover, equation ([2]) becomes 

5e-B,i = C'^ ■ {di - 5ni) , 
which, in the continuous limit leads to the Bose energies 



[l-/(ei)]dei = fceTlog 



1 + e-f^i' 



(34) 
(35) 

(36) 



From equations (fTOl) and (fTTI) follows that v = 0, so the Bose gas is also 
noninteracting, with //b = /i — ep iSj and eB = 0. 
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To test the Bose-Fermi description equivalence we invert equation (I36p 



to get 



which, if plugged into gives 



1 + e 



1 + e 



1 



(37) 



(38) 



But for a Fermi gas of constant DOS we have the identity (from equation 4, 
Ref. [8]) 

exp(/3/i) = exp(/3eF) — 1 , (39) 



so 



1 + = 1 + [exp(^eF) - 1]"^ = exp[/3(eF - //)] = exp[-/3/XB] . (40) 

Plugging (f371) into ([lOl) we obtain equation (f33l) . 

5 = exp[/3(e - = (exp{/3[eB(e) - Mb]} - 1)"' , (41) 

which proves the equivalence between the Bose and Fermi descriptions of the 
gas with constant DOS. 

The equivalence of the Bose and Fermi descriptions does not prove the 
canonical-grandcanonical ensemble equivalence. The ensemble inequivalence 
in condensed Bose systems have been extensively studied (see [El E] and 
references therein); from the EST perspective, the fluctuations of the Bose 
and Fermi energy level populations for systems of constant DOS have been 
compared in Ref. jllj . 

5.2 NON-CONSTANT DENSITY OF STATES 

Again, if we take /(e) = [exp(/?e— /?/x)+l]~"'^, then Nq = Q (no "condensation"). 
If o-(e) = Ce' with s / 0, then 



N 



and 

f^B(eB(e)) 
where again 



s + 1 



10 e/'('!-A')+l 



's+1 


1/(^+1) 


C 





g-(e)de 



l+e-''('=-A') 



1/(^+1) 



(42) 



(43) 



-/3[e(eB)-M] 



(44) 
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Using the equations (I42ll44p we rewrite the equation ([2]) as 



dee = de 



's+l 


l/{s+l) - 


C 





t cr{e) de 
l+e-/3(f-M) 



1/(^+1) 



1 + e-^^^-'^) 



M _ o-(e)de 
Jo e'9(^-A')+l 



(45) 



to obtain 



dei Cef 
1 1 + e-^(^i-'^) 



's+l 


l/(s+l) - 


C 





l+e-''("=2-M) 



1/(^+1) 



ei 



AT _ rei o-(£2) de2 
Jo e'3(^2-M)+i 



(46) 



For the two cases, e' > e and e' < e, the interaction potential u is 

1 Cef 



v{e,e' > e 
X J iV- 



dei 







de2 37 7 

1 + e^^^a-/^) 



crB(eB(ei)) 1 + e-/3(^i-/^) 



+ 







^1 Ce| dea 
1 



dei 



^1 Cei de2 

1 + e-^(^2"^) 

£1 



o-B(eB(ei)) 1 + e-/5(^i-A') 
+ 



io 1 + e-/5(^2-/^) 

s \ -1" 



7o 1 + e/^(^2-M) 
eiC^(s + 1)^ 



N + de2- 

1 + e-/^(^2-/^) 



(47) 



and 



vie,e'<e) = - r 
Jo 



dei 



Cef 



(^B(eB(ei)) 1 + e-/3(^i-^) 



(48) 



AT- / de2 ^^^^ 

1 + e'^C^a-/^) 



n -1 



. + 1) 



AT+Z de2 ^ 

^ 1 + e-/^(^2-M)_ 



-eiC=+i (s + 1)^+1 



N+ / des 37 7 

1 + e-/5(^2-^') 



dei 



Ce? 



1 _|_ g-/3(ei-/i) 

C77TT (s + l)iTT (at + /Ji dea fT(e2)/[l + e-^(^2-M)] 



s + l 



-1 



(Cef 



N - dea 



Ce2 



l_|_e/3(E2-A') 
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5.2.1 Dimensionless expressions 



To identify the relevant parameters of the gas, and for the convenience of the 
calculations, let us isolate dimensionless quantities from the equations (14211460 . 
Hereafter we shall use extensively the notations x = /3e, x' = /3e', y = and 
also their Bose correspondents xb = (3eB and = /3eg. A quantity that 
appears in almost all the expressions is the integral over Fermi distribution: 



Fi{T,C,e,n,s) = 



' Cef dei 



CikBT) 



s+1 



e^(^i^'^) + 1 
C{kBTr^'Fir{x,y,s), 



x=l3e 



x\ dxi 



(49) 



related to the integral over the hole density, 

Cef dei 



Fi'{T,C,e,^i,s) 



+ 1 



= C{kBTy+^Fir'{x,y,s) 



CikBT) 



s+1 



s + 1 



Fir{x, y, s) 



(50) 



Both, Fi and Fi' are defined for s > —1. Another quantity of interest is the 
"reduced" particle number 



N = ai±}}]:l = ( 

'-c{kBTy+^ KkBTj 



s+1 



(51) 



With the definitions (I49ll5ip . the density of states becomes 

C{kBTY [Nr-is + l)Fir{x,y,s) 



o"B(eB(e)) 



(s + 1) [Nr + {s + l)Fir'{x,y,s)] 



1/(^+1) 



CikBTY 



' is + l) 
and the Bose energy may be written as 



(^BA^riy,s),X,y,s) 



(52) 



eB = {s + 1) 



dxi xf 



1 + e-^i+2' 



• crB,r(^»-(y'S)'2;i'y'S) = {kBT)eBr{x,y,s) . 



The interaction term v{eB,e'B) (I12P has two different expressions, depending 
on the sign of ee — Cg, or equivalently, e — e': 



v{e,e > e) 



dxi xf 



s + 1 



s + l 

CikBTY Jo (1 + e-^^+y)aBr I iVr- - (s + l)Firixi,y, s 



Nr + is + l)Fir'ixi,y, s + 1) - xi [Nr + (s + l)Fir'(xi, y, s)] 



s/{s+l) 



(53) 
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and 



v{e,e' < e) 



^ + 1 
CikBTy 



dxi xf 



:i+e-^i+S/)(TBr 
1 



s + 1 



Nr - {s + l)Fir{xi , y, s) 



Nr + {s + l)Fir'{xi,y, s) - xi [Nr + {s + l)Fir'{xi,y, s)] 



+ 



s + l 

CikBTy 



dxi 



[Nr + (s + l)Fir'(xi,y,s)]-+i 



1 + e-^^+y 



Nr - {s + l)Fir{xi,y,s) 



^■(54) 



5.2.2 Low energy expressions for s ^ 



All the formulas in Section 15.2.11 are very difficult to calculate numerically. 
Therefore, low energy asymptotic expressions can be very useful in calcula- 
tions. To find the low energy behavior of the functions involved, we calculate 
the Taylor expansions of the dimensionless quantities around e = (x = 0). 
We start with the integrals over Fermi distributions. 



1 



X 



1 



^ xf dxi 

e'^i-y + 1 ^ e-y + 1 s + l~ (e-y + 1) • (e!^ + 1) s + 2 



(55) 



xf dxi 



1 



s+l 



+ 



Fir{x, y, s) = 
and 

Fir'{x, y, s) 



for X ^ 1. Using Eqs. (j55p and (j56p we get asymptotic expressions for all the 
other variables. The Bose DOS is 



e?^-^i + 1 ey + ls + l (e-y + 1) (e?^ + 1) s + 2 
+1 

Fir{x, y, s) , (56) 



s + l 



'7B,r(a;,y,s) 







A^^ - 


x^+l 1 








e-y+1 






1 + 


1 




X 


s+l ■ 


Nr{ey+1) Ari/(s+i) 



(57) 



e-y + 1 



1 + 



-y 



s + l 

where 6[x) is the Heaviside step function. At x = xb = 0, (TB,r = Nr^^^~^^\ 
so it attains a finite value, even if a is either zero or infinite. Nr^^^^^^ is the 
average density of the Fermi gas, from zero to ep. If s > 0, aB,r is increasing 
with X, while for s < 0, aB,r is decreasing. This is a general result, since already 
from equation ([3]) one can observe that o"b is monotonic for monotonic a. 
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Bose energies are: 



ee « {s + l)kBT r - 
Jo i 



dxi xf 


AT _|_ ^1 
^^r -t- eM+l 


1/(^+1) 

- Xl 


1 + e-^i+f 




A^^ - 


] 








e~y+i 





kBTNr 



-s/{s+l) 



+ 



e^ + 1 
2 iV^(e-2^ + r 



x"+^ - 9(s) ■ x' 



+2 



S + 1 



s + 2 \iv:y(^+i) e-S' + 1 



• 1 + 



(58) 



In the lowest order, for both, s < and s > 0, xb have the same expression: 

6B N-'/^'^'^ 



XB 



k^T ev + 1 



(59) 



By ehminating x from Eqs. ()59p and ()57p we obtain the asymptotic expression 
of crB(eB), 



,,(x,y,s) « iV;/(^+i)+0(s)-iV/^(eJ^ + l)^/(^+i)x^/('+') 



-s) -iVZ+^e^ 



1 + 



s + 1 



XB . 



(60) 



So cJB,r decreases hnearly with xb for s < 0, while for s > it has an infinite 
positive slope at xb = 0. 

The expressions for the interaction potential are more complicated. Let 
us start with x' > x. After some calculations we arrive to the following results: 



t;(x, x' > x) 



s + 1 Nr 



C{kBTy ev + 1 
for s > 



+ 



s + 1 



s + 1 ' s + 2 \ 1 + e-?^ j^TTT 



and 



v{x, x' > x) 



-2s-l 
Nr 



C{kBTy eV + 1 



sx^+i + 



2s + 1 



1 + 



2(s + l) I s + 1 



X 



2s+2 



Nr{e~y + 1) 



for s < . 



Obviously, v{x,x' > x) — > 0, as x ^ 0, since eB = is the ground state for 
any configuration. 
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Bosonic "quasiparticle" energy in 1 D Bosonic "quasiparticle" energy in 3D 




2468 2468 10 



Figure 4: The Bose quasiparticle energy, ee, vs. e, for (a) ID {s = —1/2) and 
(b) 3D (s = 1/2) systems, and the Bose DOS, (TB(eB) for the same systems. 
In ah the cases we take ep = lOk^TQ, where Tq is a scahng temperature. The 
four curves in each pfot correspond to four values of T: O.ITq (solid line), Tq 
(dotted line), 2To (dashed line), and ATq (dash-dot line). 



For X > x', 
v{x, x' < x) - 



s + l 



Nr 



-2s-l 

3 + 1 



C{kBTyNr ey + 1 



{x - x') + 



J\2 



X' 



s+l 



2{l + e-y) s + l 



l\s+l 



, for s > 0, 



(61) 



and 



v{x, x' < x) 



C{kjiTY{ey + 1 



{x'y+^ - (X - X') -{s + l) 



2{s + l)Nr{l + e-y) 
for s < 0. 



s + l- se-y) - s{x'y+^ - ^^i^ — ^T!^ 



2{s + I) Nr{e-y + I) 



(62) 



In both equations (|61|) and (|62|) only the two lowest orders are correct, whichever 
these are, depending on s. Again, note that t>(x, x' < x) ^ as x — > 0, for 
any s. 

In Figs, m (a) and (b) it is shown the dependence of ee on e for ID 
and 3D gases and four different temperatures. The dependence of ctb on ee is 
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fee, 
bp 





(a) 











e/ikBT) 




e/{kBT) 



Figure 5: For a ID Fermi gas (s = —1/2): (a) The population of the bosonic 
energy levels, 6F[eB(e)] and 6B[eB(e)], for y = /^/(/cbT) = 10 (1) and y = 
^/{k-QT) = 20 (2). For both values of y, the curves corresponding to ftp and 
6b are indistinguishable, (b) The bosonic effective quasiparticle energy, eB(e), 
for the same values of y as in (a). 



shown in figured] (c) and (d) for the same gases and temperatures. The four 
curves in each of the plots correspond to four different temperatures, namely 
O.ITq, To, 2To, and ATq, where Tq is a scaling temperature and ep is chosen so 
that ep = 10A;bTo. 



5.3 THE EQUIVALENCE OF THE BOSE AND FERMI DE- 
SCRIPTIONS 

Having now the bosonic picture of our Fermi system, built based on the Fermi 
distribution, we can analyze the bosonic ensemble as in Section [H Using the 
results of Section [5.21 we calculate b by both formulas (|24p and (j33p and com- 
pare the results. To avoid confusion, we denote by 6p and 6b the populations 
given by equations (p^ and (j33|) . respectively. Since the two systems are 
thermodynamically equivalent, 6p and 6b should be identical in the thermo- 
dynamic limit. This identity appears to be difficult to prove analytically, so 
we calculated numerically eB(e) and with this, 6b, which we then compared 
with 6p. As one can see in figure [5l 6p[eB(e)] and 6B[eB(e)] are equal within the 
numerical accuracy, for any ^/{k^T). (Note: we plotted the numerical results 
for s = —1/2 because for this value of s the results are much more accurate 
and the multiple integrals involved converge for a wider range of chemical 
potentials and single particle energy levels.) 
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6 CONCLUSIONS 



In this paper it is introduced the basic idea and results of what have been 
caUed exclusion statistics transformation (EST). The EST from Fermi to Bose 
systems is done in the following way: the fermions are pictured as piles of 
particles that occupy completely intervals along the single particle energy axis 
(no empty single particle states in each of the intervals) , which are then rotated 
into horizontal position by associating to all of them a single "Bose" energy 
(see figure [1]). Using EST we transformed an ideal Fermi gas into a Bose 
gas with the same excitation spectrum. Having the same excitations, the two 
gases are thermodynamically equivalent by construction [6], i.e. if the Fermi 
gas has entropy Sp, internal energy Up and particle number N, then the Bose 
gas has the same entroy at any temperature as the Fermi gas, Sb = Sp. The 
internal energy of the Bose gas is Ub = Up — C/g.s, where C/g.s. is the energy of 
the Fermi system at zero temperature. 

Having constructed in this way the Bose gas, one can calculate all the 
thermodynamic quantities of the Bose and Fermi gases independently, by max- 
imizing the partition functions with respect to the quasiparticle levels popula- 
tions. Based on the thermodynamic equivalence, all the canonical properties 
of the two gases should be identical and, moreover, applying EST, the Fermi 
distribution should transform into the Bose distribution and vice-versa. We 
checked this identity numerically, by calculating the Bose and Fermi popula- 
tions. The populations transform indeed into each-other by EST within the 
numerical accuracy. So the thermodynamic equivalence is not applicable only 
to a very special class of systems, namely ideal systems with the same, con- 
stant, DOS, but should be regarded as a very general concept. We have now a 
method to transform gases of a given exclusion statistics into equivalent gases 
of some other exclusion statistics. 

Although this is not discussed here, at low temperatures the Bose system 
condenses and the Bose condensate is related by EST to the particles in the 
Fermi system that stay on the lowest energy levels and occupy completely an 
energy interval. By analogy to the Bose system, these particles form the so 
called Fermi condensate [9]. Using this analogy, the Fermi system may be 
described as consisting of a condensate of, say A'^o particles, and with N — Nq 
particles above the condensate. The condensate forms in interacting systems 
Ref. [8l Ho]), but also in ideal systems [9j. At low temperatures, when the 
condensate lies close the the Fermi energy, the single particle spectrum of the 
"thermally active" N — Nq particles may be approximated as constant and the 
EST becomes very simple (in any number of dimensions), since it reduces to 
the EST in systems of constant DOS 



It is well known from Bose systems that the properties of the condensate 
are strongly influenced by the dimensionality of the system (see for example 
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[15j). but the collective properties of the condensate, which certainly deserve 
investigation, were omitted in this paper. 
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